	
	*****************************************************************
	* Taxing the Poor Twice: Bandwidth and the Value of Consumption * 
	*****************************************************************
	
	* Figure S1 Consumption Choices 
	
	
	
	* 1. Importing Merged Data 

	use "$dir_data/bandwidth_long_final.dta" , clear
	
	* setting scheme for figures
	cd "${bandwidth}"
	set scheme eop

	keep pid treat treated choice_song choice_video choice_food choice_activity 
	collapse (firstnm) treat treated choice_song choice_video choice_food choice_activity, by (pid)


	
	* 2. Figure 

	* Choice of song 

mlogit choice_song i.treated
test 1.treated

	// p-value=0.1864

	* graph margins 
margins treated, atmeans predict(outcome(1))
marginsplot, name(song1_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Song 1", size(large))

margins treated, atmeans predict(outcome(2))
marginsplot, name(song2_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Song 2", size(large))

margins treated, atmeans predict(outcome(3))
marginsplot, name(song3_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Song 3", size(large))

margins treated, atmeans predict(outcome(4))
marginsplot, name(song4_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Song 4", size(large))

graph combine song1_p song2_p song3_p song4_p, ycommon imargin(6 6 6 6) title("Song Choice Probabilities by Treatment Status")

graph export "$dir_figures/figureS1_song_choice_probabilites.eps", replace 


	* Choice of video

mlogit choice_video i.treated
test 1.treated

	// p-value=0.1920

	* graph margins 
margins treated, atmeans predict(outcome(1))
marginsplot, name(video1_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Video 1", size(large))

margins treated, atmeans predict(outcome(2))
marginsplot, name(video2_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Video 2", size(large))

margins treated, atmeans predict(outcome(3))
marginsplot, name(video3_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Video 3", size(large))

margins treated, atmeans predict(outcome(4))
marginsplot, name(video4_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Video 4", size(large))

graph combine video1_p video2_p video3_p video4_p, ycommon imargin(6 6 6 6) title("Video Choice Probabilities by Treatment Status")

graph export "$dir_figures/figureS1_video_choice_probabilites.eps", replace 

	* Choice of food  

mlogit choice_food i.treated
test 1.treated

	// p-value=0.0427

	* graph margins 
margins treated, atmeans predict(outcome(1))
marginsplot, name(apple_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Apple Slices", size(large))

margins treated, atmeans predict(outcome(2))
marginsplot, name(chocolate_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Chocolate", size(large))


margins treated, atmeans predict(outcome(3))
marginsplot, name(samosas_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Samosas", size(large))

graph combine apple_p chocolate_p samosas_p, ycommon imargin(6 6 6 6) title("Food Choice Probabilities by Treatment Status")

graph export "$dir_figures/figureS1_food_choice_probabilites.eps", replace 

	* Choice of activity 

mlogit choice_activity i.treated
test 1.treated

	// p-value=0.9131

	* graph margins 
margins treated, atmeans predict(outcome(1))
marginsplot, name(act1_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Activity 1", size(large))

margins treated, atmeans predict(outcome(2))
marginsplot, name(act2_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Activity 2", size(large))

margins treated, atmeans predict(outcome(3))
marginsplot, name(act3_p) xlabel(0 "Control" 1 "Treated", labsize(large)) title("") xtitle("") ytitle("Prob. of Activity 3", size(large))

graph combine act1_p act2_p act3_p, ycommon imargin(6 6 6 6) title("Activity Choice Probabilities by Treatment Status")

graph export "$dir_figures/figureS1_activity_choice_probabilites.eps", replace 


